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ABSTRACT 

Spectroscopic, photometric and dynamical data of the inner 3 kpc part of the starburst galaxy M 82 
are analyzed in order to investigate the star formation history of its disk. The long-slit spectra along 
the major axis are dominated by Balmer absorption lines in the region outside the nuclear starburst 
all the way up to w 3.5 scalelengths (/xg = 22 magarcsec -2 ). Single Stellar Population (SSP) spectra 
of age 0.4-1.0 Gyr match well the observed spectra in the 1-3 kpc zone, with a mean age of the 
stellar population marginally higher in the outer parts. The mass in these populations, along with 
that in the gas component, make up for the inferred dynamical mass in the same annular zone for 
a Kroupa initial mass function, with a low mass cut-off mi = 0.4 M© . The observed ratio of the 
abundances of a elements with respect to Fe, is also consistent with the idea that almost all the 
stars in M82 disk formed in a burst of short duration (0.3 Gyr) around 0.8 Gyr ago. We find that 
the optical/near infrared colors and their gradients in the disk are determined by the reddening with 
visual extinction exceeding 1 mag even in the outer parts of the disk, where there is apparently no 
current star formation. The disk-wide starburst activity was most likely triggered by the interaction 
of M82 with its massive neighbor M81 around 1 Gyr ago. The properties of the disk of M82 very 
much resemble the properties of the disks of luminous compact blue galaxies seen at 0.2-1.0 redshift. 
Subject headings: galaxies: individual (M 82) — galaxies: evolution — galaxies: interactions 



1. INTRODUCTION 

M82 is a nearby edge-on galaxy (D = 3.63 Mpc, 
image scale 17.6 pcarcsec -1 ; Freedman et al. 1994), 
which harbors a nuclear starburst, the prototype of the 
starb urst phenomenon l)Rieke et al. 1 119801 iRieke et"aTl 
1993). The optical appearance of M82 is dominated 
by brig ht star-forming knots interspe rsed by dusty fil- 
aments Ip'Connell fc Mangano lll978h . This appearance 
has led to an Irr II morphological classification of M 82 
(Hol mberg 1119501) . A near- infrared (NIR) b ar of ~1 kpc 
length was discovered bv iTelesco et al. I l)1991|) . Re- 
cently, [Maj^j^£amiic^&^unjy||2005) have discovered 
two symmetric spiral arms, suggesting a morphological 
classification SBc. Embedded inside the 500 pc star- 
burst region is a small bulge, with an upper ma ss limit 
of 3 x 1O 7 M llOaffnev. Lester fc Telescol Il99l . The 
if-band scalelen gth o f the disk is as small as 860 pc 
HTchikawa et al. BlflflfiTL 

M 82 is a low mass galaxy (~ 10 10 M Q ), with a strong 
central concentration of mass and har dly any evidence 
for the dark matter halo l|Sofue 1119981) . This' mass is m 
the lower mass range of late-type galaxies. However, the 
observed luminosity, mean surface brightness and metal- 
licity are higher than expected for a late-type galaxy. In 
addition, the percentage of mass in the gaseous compo- 
nent is as much a s ~ 50%, which is abnor mally high for 
a normal galaxy fYoung fe Scoville"lll984[) . Several au- 
thors have attributed these abnormalities to some phe- 
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nomenon associated wi th its interact ion with other mem- 
bers of the M 81 group. iSofue I (^98) explained the small 
size of the disk and the falling rotation curve with a 
model wherein M82 lost its outer disk and halo during 
a closeby interaction with M 81 and NGC 3 77, w hich 
started around 1 Gyr ago. lYun. Ho. fc Lo I l)1994fi re- 
produced the observed structure of the HI streamers 
using a numerical simulation, where M 82 is interacting 
with M 81 and NGC 3 077. The close st approach occurred 
around 0.3 Gyr ago. lElvius 1 1)1972|) presented a picture 
wherein M 82 captured a gas cloud during its passage in 
the M81 group, in order to explain the high gas frac- 
tion. Are the observed properties in the disk of M82, 
especially, its high surface brightness, high metallicity, 
and low mass-to-light ratios, consistent with these sce- 
narios? In this paper, we investigate the star formation 
history (SFH) that is consistent with the observed prop- 
erties of the stellar disk, and compare the resultant SFH 
with phenomena that may be related to the interaction. 

In Section 2, we describe the spectroscopic and imag- 
ing observations used in this work. Section 3 deals with 
a discussion of the quantities that can be used as con- 
straints to the age of the stellar populations in the disk 
of M82. In Section 4, we discuss the star formation and 
chemical evolutionary histories that are consistent with 
our data. The evolutionary status of M 82 is discussed in 
Section 5. A summary of the main results of this work 
is given in Section 6. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Spectroscopic Observations 

Spectroscopic observations of M82 were carried out 
with the purpose of registering the absorption-line domi- 
nated stellar spectra at a variety of radial distances from 
the center. Long-slit (3' in length) spectra were taken at 
three positions along the major axis of the galaxy (P.A. 
of the slit = 62°) using the Boiler & Chivens spectro- 
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graph at the Guillermo Haro Astrophysical Observatory, 
Cananea, Mexico. The observations were carried out on 
three nights, one in 1999 February and the other two in 
2003 November. The starburst nucleus was placed at the 
center, extreme southwest and extreme northeast of the 
slit, respectively on the three nights. With the three slit 
positions, spectra up to a major axis distance of 3' were 
registered on either side of the nucleus. At this radial 
distance, £?-band surface brightness levels correspond to 
22 magarcsec -2 . Two spectra, each of 30 minute du- 
ration, were taken for the bright central part, whereas 
three 30-minute spectra, were taken for the other two 
slit positions. A slit width of 1.6" was employed. A grat- 
ing of 150 lines/mm was used which covers a wavelength 
rang eof «3500-6800A at a spectral resolution of around 
10 A. The spectral and spatial samplings were 3.2 A/pix 
and 0.46"/pix respectively. Object was observed as close 
to the meridian as possible, with maximum airmass of 
« 1.30. The instrumental response was calibrated by the 
observation of standard stars HR1544, HR5501, BD+40 
4032, Feige 15 and Feige 34. 

Each frame was bias-corrected and divided by a nor- 
malized flat-field using the tasks in the IRAF package. 
Wavelength-calibrated frames of the same slit position 
were averaged, in the process removing cosmic ray events. 
The disk of M 82 occupies the entire slit-length and hence 
sky spectrum could not be extracted from the same spec- 
trum as the object. Instead, sky spectra were extracted 
from the spectra of galaxies with a compact nucleus ob- 
served immediately before and after the M82 observa- 
tions. 



and southwest directions separately. Beyond a radius 
of 1 kpc on either side of the nucleus, bright knots are 
absent, and the spectra represent that of the underlying 
smooth disk. We defined two spectra on either side of the 
nucleus corresponding to radial zones between 1-2 kpc, 
and 2-2.7 kpc. The positions along the major axis of 
the extracted spectra are identified in Figure 1 by verti- 
cal hatches, which are labeled by their mean distances in 
arcsecs from the starburst nucleus. In this figure, we also 
show the major axis intensity cuts in the optical and NIR 
bands (solid lines). The exponential disk profiles (ob- 
tained by fitting the azimuthally averaged intensity pro- 
files) are also plotted in this figure (dashed lines). The 
bottom two panels show the major axis color profiles. 
The dashed vertical lines mark the position of brightest 
region of the recently-discovered spiral arms, where the 
colors are marginally bluer. It is easily noticeable that 
colors become systematically redder towards the center 
of the galaxy. 

The most critical part of the spectrum for population 
synthesis analysis is the blue part of the spectrum, be- 
cause of the presence of a variety of spectral features 
sensitive to young, intermediate-age and old popula- 
tions. This part of the extracted spectra are plotted 
in Figure 2. The signal-to-noise ratio of these spectra 
are good enough to measure spectral indices sensitive 
to age and metallicity. Some of the important absorp- 
tion lines are marked in the figure. In general, all spec- 



2.2. NIR and Optical Images 

NIR images in the J, H and K bands used in the 
present work were taken with the CAnanea Near In- 
frared CAmera (CANICA) 5 available at the 2.1-m tele- 
scope of the Observatorio Astroffsico Guillermo Haro 
in Cananea, Sonora, Mexico. Details of these images 
and the analysis procedure ado pted were described in 
iMavva. Carrasco &: Luna 1 1)2005|) . Optical images in the 
U, B, V and R bands are based on the data archived at 
the NASA Extragalactic Database. These images were 
taken at the 40-inch telescope of the Mt. Laguna O bser- 
vatory and were described bv lMarcum et alJ lj2001|) . 

3. CONSTRAINTS ON THE DISK PARAMETERS FROM 
OBSERVATIONAL DATA 

In this section, we combine our spectroscopic and pho- 
tometric data with published dynamical data to put 
strong constraints on critical parameters of the stellar 
disk of M82. 

3.1. Spectroscopic ages of the stellar populations 

Two-dimensional spectra of M 82 were closely exam- 
ined to isolate all regions with noticeably distinct sur- 
face brightness values. On either side of the starburst 
nucleus, there are bright continuum emitting knots up 
to around 60" (1 kpc) distance (e.g. region B on the 
northeast and region G on the southwest as defined by 
O'Connell & Mangano 1978). Spectra representative of 
these continuum knots were extracted on the northeast 



see http://www.inaoep.mx/~astrofi/cananea/canica/ 
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Fig. 1. — One-dimensional intensity and color profiles of M82. 
In the top panel, we show the intensity cuts along the major axis 
in K,J,R, and B bands (solid lines), along with the exponential 
function that fits best the azimuthally averaged intensity profiles 
in the corresponding bands. The positions of the extracted spectra 
are shown by vertical grids, with the associated numbers denoting 
the mean distance in arcsecs of extracted regions from the nucleus. 
The NE and SW directions are also shown. The bottom two pan- 
els depict the major axis color profiles. The vertical dotted lines 
denote the intersection of the spiral arm with the major axis. 
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Fig. 2. — Major-axis spectra of the stellar disk of M82 at five 
selected zones. In (a), the top two spectra correspond to the north- 
east (NE), and the bottom three correspond to the southwest (SW) 
part of the disk. The numbers appearing after the NE or SW cor- 
respond to the radial distance from the nucleus in arcsecs of the 
extracted spectra. Prominent age-sensitive absorption lines are 
identified. Note that the almost feature-less part of the spectra 
between H7 and H/3 is not shown, (b) The 60" and 110" spectra 
from either side are superposed (line types and colors the same as 
in panel a). Note that the strengths of all the absorption features 
in the four spectra are nearly identical, despite the NE spectra 
being systematically bluer than the SW spectra. 



Fig. 3. — Observed spectra (dotted blue line) compared with 
model SSP spectra (solid red line), (a) Six models are superposed 
on the observed SW60" spectrum. The models correspond to two 
metallicities (Z=0.050, 0.020), and three ages (0.3-0.6 Gyr). The 
best-match is obtained for 0.5 Gyr at Z=0.020. (b) The best- 
match model SSP superposed on the rest of the inner disk spectra. 
All these three spectra are consistent with an SSP of 0.5 Gyr at 
Z=0.020, with the only difference between the spectra being the 
amount of reddening, (c) Three models are superposed on the SW 
160" spectrum. This spectrum differs from the inner disk spectra 
in having the age-sensitive Ca II K profile almost as deep as the 
Ca II H+He profile, an indication of relatively older population. 
The model that best matches the important observed features is 
at Z=0.020 for 0.9 Gyr. 



tra are dominated by Balmer absorption lines, which 
is consistent with the "A-F" spectral type inferred by 
lO'Connell fc Mangano I l)1978|K Relative depths of the 
spectral lines in the four spectra corresponding to the 
regions in the inner disk are identical. This is illustrated 
in panel (b) of the figure where these four spectra are su- 
perposed. The only difference is in the continuum shape, 
with the NE spectra systematically bluer than the SW 
spectra. In comparison, the spectrum of the outer disk 
(SW110") departs from the inner disk spectra in having 
deeper Call K line as compared to the Call H+He and 
other Balmer lines. 

In the remaining part of this section, we use the ob- 
served spectra to infer the age of the stellar populations 
that give rise to these features, using model spectra for 
Single Stellar Populations (SSPs), properly smoothed to 
the observed resolution. The model SSPs have a Kroupa 
like initial mass functio n (IMF) — with masses be- 
tween 0.15 and 120 ilKroupa et al1ll993l). They are 
based on the Padova l|Bressan et al.lll993t |F^ggJJo^t^L 
1994a|3) evolutionary tracks and on the l.lacobv et al. 
( 1984) library of stell ar spectra. More d etails of the 
model can be found in lBressan et al.l l)1998[) . It is impor- 
tant to stress here that while we compare observed spec- 
tra with SSP models computed with the usual Kroupa 
IMF, in Sec. 4 we will discuss also a top-heavy modifica- 



tion of the IMF. The use of such a top-heavy IMF does 
not change our estimate of the age of the disk of M82. 

Reddened model spectra are superposed on the ob- 
served spectrum for the SW60" zone in Figur e 3a, where 
we have made use of ICardelli et al.l (|1089f) reddening 
curve. The metallicity, age and the visual extinction of 
each model is indicated below each displayed spectrum. 
The intensity ratios of different features in the observed 
spectrum of the inner disk can be reproduced only for 
ages between 0.3 and 0.6 Gyr, with the exact age de- 
pending on the choice of metallicity. For ages younger 
than 0.3 Gyr, the Ca II K line does not appear in the 
SSP, whereas for older ages it is much stronger than the 
observed value. The model that best matches the ob- 
served spectrum corresponds to an age of 0.5 Gyr and 
solar metallicity. In Figure 3b, we illustrate that all the 
spectra in the inner disk be reproduced at solar metallic- 
ity for an age of 0.5 Gyr, with extinction being different 
at different positions of the disk. Clearly, the southwest- 
ern side of the disk is more reddened as compared to 
the northeastern side. Derived ages would be younger 
(older) by around 0.1 Gyr if the metallicities are higher 
(lower) than solar values. More importantly, our spectra 
of inner disk regions are not consistent with mean ages 
greater than 0.7 Gyr. 
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The outer disk spectrum suggests a population that 
is marginally more evolved as compared to that of the 
inner disk. Nevertheless, the 4000 A break, which is a 
signature of populations a few gigayears old, is absent. 
We estimate an age of 0.9 ±0.1 Gyr for the outer disk. 
In principle, acceptable fits can be obtained by adding 
an old population of several gigayears, typical of galactic 
disks, to a population of 0.5 Gyr old. However, single 
SSPs give a better match to the observed spectrum than 
those obtained by adding an underlying old population. 

As noted earlier, the NE60" spectrum corresponds to 
the the fossil starburst region B. The Hubble Space Tele- 
scope (HST) images of region B have revealed that this 
region c onsists of more than hundred super star clusters 
(SSCs). IdeGriis. O'Connell fc Gallagher! lEffll have es- 
timated a mean age of these SSCs to be 0.65 ± 0.25 Gyr, 
with none of the clusters younger than 0.30 Gyr. Typ- 
ical visual extinction values were found to be less than 
1 mag. These ages were estimated based on the optical- 
NIR (BVIJH) spectral energy distributions of individual 
clusters, keeping the reddening E(B-V) as a free parame- 
ter. Our spectroscopically estimated age of 0.5 ±0.1 Gyr, 
and A v — 0.5 mag are very much consistent with 
these photometrically estimated values. However, a re- 
analysis of the s ame HST dataset has yielded older ages 
1.10 ± 0.25 Gyr <|de Griis. Bastian fc Lamersll2003D . but 
with mean A v < 0.2 mag. It is well-known that both age- 
ing and reddening have the same effect on the observable 
photometric quantities, and hence the age and extinction 
cannot be derived completely independent of each other. 
On the other hand, in the spectroscopic method that we 
have adopted, the relative strengths of absorption lines 
fix the age, whereas the slope of the continuum fixes the 
extinction, thus the age, and extinction were determined 
independent of each other. We re-iterate that the contin- 
uum of our spectra are consistent with both the solutions 
(low age with high extinction or high age with low ex- 
tinction), but the absorption- line ratios clearly rule out 
ages exceeding 0.7 Gyr for region B. 

The similarity of the inner disk spectra, irrespective 
of whether they belong to compact knots of region B or 
diffuse parts of the disk, carries important clues on the 
star formation scenario in the disk of M82. In the next 
section, we use the mass-to-light ratio in the if -band to 
establish firm constraints on the mass of the underlying 
old disk population, and also discuss a scenario of star 
formation in the disk of M82. 

3.2. Absolute magnitude and colors 

Azimuthally averaged intensity profiles of M 82 in the 
NIR and the optical bands follow an exponential law ex- 
pected for stellar disks of galaxies (see Figure 1). We 
generated the 2-dimensional exponential disk of M82 
in the if-band based on the observed scalelength of 
47" (825 pc). We then extracted the disk magnitude in 
the annular zone between 70" and 170", which was found 
to be 5.84 mag. At the distance of M 82, this corresponds 
to an absolute magnitude of —21.96, without correction 
for extinction. 

Optical and NIR colors are useful in determining the 
evolutionary status of the stellar systems, if the redden- 
ing in various colors can be reliably estimated. The use 
of multiband images allows us to obtain profiles of az- 
imuthal and radial distribution of colors. In Figure 4, we 
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Fig. 4. — Azimuthal color profiles at radius=110, 120 and 130" 
from the center. The position of the spiral arms is marked by 
the vertical lines. The arms, in general, are bluer than the disk. 
The disk colors in the azimuthal zone between 65-245° is affected 
by the foreground dust and the superwind cone, and hence typical 
disk colors and their errors are derived averaging the colors between 
P. A. 320° and 360°. 



display the azimuthal distribution of colors in the radial 
zone between 90" and 130". In this radial zone, the arms 
are the strongest features, and hence these profiles can 
be used to investigate whether the results obtained from 
spiral-arm dominated major axis profile is applicable to 
the rest of the disk. In general, colors are the bluest 
along the major axis (dotted lines at P.A.s 65 and 245°). 
However, the color difference between the disk and the 
spiral arms at a given radius is small compared to the 
overall radial gradient in colors, displayed in Figure 1. 
Hence it seems that the results obtained from major axis 
spectra are representative of the entire disk. 

3.3. The Mass and mass-to-light ratio 

The rotation curve of M 82 has been derived in several 
studies using both stellar and gaseous tracers up to a 
radius of 170", with th e intention of mappin g the mass 
distrib ution in its disk. Ma vall I i|1960l) and lGoetz et alJ 
(1990) obtained a rotation c urve using op tical emission 
and absorption lines, whereas fSofue 1 1)1998|) used the CO 
and HI lines. In all these studies, the rotation curve has 
been modeled to obtain the radial distribution of mass. 
iMavall I l|196(y estimated a mass of 4 x 10 9 M Q in the 
70-170" an nulus, whereaslGoetz et alJ (|1990ft masses are 
50% lower. iSofue I l|1998fl found that most of the mass 
is concentrated within the central 1 kpc radius, with the 
mass distribution outside this radius consistent with an 
exponential mass surface density profile. Significantly, 
there is no evidence for a dark matter halo, and hence 
the entire mass outside the central bar can be associ- 
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ated with the stars and gas in the disk. We obtained 
the mass-to-light ratio 6 of the disk outside a radius of 
1 kpc by combining the disk mass and the -ftf-band lumi- 
nosity in the same annular zone. The resulting value is 
0.30 ± 0.05 Mq/Lkq if we use mass estimat es of iM avall 
(1960) and 0.15 for the mass estimates of iGoetz et al.l 
(1990). At least 15% of the disk mass in the annular 
zone ou tside 1 kpc is in the form of atomic and molecu- 
lar gas (jYoung fc Scoville 111 984J). 



Table 1. Observed and model properties of the disk of M82 



Quantity 


units 


Observed 


Model A 


Model B 


A^disk 
M(gas) 
M(disk) 

M K (disk) 

M/L K 

[O/Fe] 


1O 9 M 
mag 

M /l Kq 


< 4.0 
> 0.15 
-21.96 ±0.05 
< 0.30 ± 0.05 
0.35 ± 0.20 


4.0 
0.15 
-21.42 
0.53 
0.11 


4.0 
0.15 
-21.95 
0.32 
0.21 



3.4. Metallicity of the nebular gas and stars 

Metallic abundances of gas and stars in the star- 
burst region of M 82 have been analyzed by Origlia ct al. 
( 2004) . The cold phase gas abundances were derived us- 
ing the optical and NIR nebular lines associated with the 
giant HII regions, whereas the stellar abundances were 
obtained using the NIR absorption lines that originate 
in the atmospheres of cool red supergiants. Since these 
stars live only for a few tens of million years, both the 
stellar and gaseous abundances are expected to represent 
the values present in the disk before the onset of nuclear 
starburst. 

In Table 1, we summarize all the observed quantities 
for the disk, which we will use as a constraint for a de- 
tailed modeling of the star formation and chemical his- 
tory of the disk of M 82. 

4. STAR FORMATION AND CHEMICAL HISTORY OF THE DISK 

OF m82 

In order to reproduce the observed quantities (metallic- 
ity enhancement, M/Lk, age and colors), we ran a num- 
ber of chemical evolutionary models from the GRASIL 
WEB interface GALSYNTH 1 . In these models, gas 
phase metallicities of different elements are calculated 
over the Hubble time for a given SFH of the galaxy. 
At each epoch, the star formation rate (SFR) is cal- 
culated using a Schmidt law, with the exponent value 
of 1.0. The standar d Kroupa IMF parameters are used 
l)Kroupa et alJ f 19931. The current M/Lk value is calcu- 
lated as the ratio of the mass locked up in stars (both 
living as well as their remnants) to the luminosity of all 
the living stars. It is important to note that in order to 
reproduce the observed M/Lk values, it is necessary to 
correct the model values for the mass in gas (> 15%), 
and the extinction in the X-band. 

Among the observed quantities, the metallicity en- 
hancement critically depends on the slope and the upper 
cut-off mass of the IMF and the burst duration, whereas 
the lower cut-off mass of the IMF and the age of the burst 
determine the M/Lk ratios and colors. The observed set 
of parameters could be produced only by those models 
where almost all the observed stellar mass of the disk was 
formed by a burst lasting for only a few hundreds of mil- 
lion years, the burst having completely stopped around 
0.5 Gyr. We ran a set of models for a wide range of input 
parameters. We discuss two of these models below. 

6 Throughout this work, mass-to-light ratio is defined using the 
monochromatic i\-band luminosity calibrated against the absolute 
magnitude of sun Mkq = 3.33 mag. Instead, if we use the _ft"-band 
luminosities integrated over a standard if-band response curve of 
band width 0.35^tm, and express mass-to-light ratio in Mq /~Lq , 
then the quoted values have to be multiplied by 50. 

7 http: / /web. pd.astro.it/galsynth/index.php 



The run of the SFR of the best-fit model as a function 
of the age of the living stars is shown in Figure 5a. In Fig- 
ure 5b, we show the metallicity enhancements for various 
metals for the standard Kroupa IMF (open circles de- 
noted model A), along with the observed enhancements 
in the cold gas and stars. The corresponding values for 
the M/Lk and colors are shown in the bottom panel 
of the figure. With the Kroupa IMF the model pro- 
duces metallic enhancements which are systematically 
outside the observed values even after taking into ac- 
count the errors in the observed metallicities. The model 
also produces M/ L K values which are around 30% higher 
than the highest values allowed by the observations. The 
model colors, after reddening them with the inferred ex- 
tinction values from the spectra, are consistent with the 
observations. 



The only way we could reproduce the observed high en- 
hancements of a elements and a low M/Lk is by adopt- 
ing a top-heavy IMF (model B) . In this model, the slope 
of the Kroupa IMF for the massive stars is reduced from 
1.7 to 1.4 (i.e. closer to the Salpeter slope), and the lower 
mass cut-off raised from 0.15 to 0.4M Q . The resulting 
values are shown in Figure 5 by square symbols. This 
model successfully reproduces the observed trend in the 
metallic enhancements for all the elements except Car- 
bon. The M/Lk values are also very well reproduced. 
In fact, the agreements between the model and the ob- 
servational values would be even better if we use a flatter 
IMF. These trends for the disk IMF are very similar t o 
that found for the nuclear IMF bv IRieke et all l)1993|) . 
The shallower IMF would imply a harder ultraviolet ra- 
diation field, contrary to what is observed fo r a sample 
of 23 starburst nuclei l|Rigbv & Rieke 1 12004^1 . This ap- 
parent contradiction could be solved if the most massive 
stars spend most of their lifetimes embed ded in dense 
highly extinguished regions, as suggested bv lSilva et al. I 
( 1998). The model quantities are compared with the ob- 
served quantities in Table 1. 

In summary, models with more or less constant SFR 
over the Hubble time cannot reproduce neither the ob- 
served low mean age of the disk of M 82, nor the observed 
chemical pattern and M/Lk values. On the other hand, 
observations can be only reproduced if more than 90% of 
the stellar mass inside the central 3 kpc is younger than a 
gigayear. Hence, around 1 Gyr ago, the disk experienced 
an intense starburst. Furthermore, we can reproduce 
both the low M/Lk value and the enhancement of the a 
elements only by a suitable modification of the Kroupa 
IMF. To reproduce the observed M/Lk, half a gigayear 
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Fig. 5. — (a) Star formation History of the disk of M82 that 
reproduces best the observed properties, (b) Observed stellar and 
nebular metallicities compared with the gas metallicities of our 
best-fit model. Model A and B correspond to standard and top- 
heavy Kroupa IMFs, respectively, (c-d) B — V and V — K colors 
are plotted against the K-band mass-to-light ratio of the disk. The 
error bars on the colors are an indication of the observed color gra- 
dient between 80-160" part of the disk. Open circles and squares 
correspond to Model A and B, respectively, without taking into 
account the effect of the gas mass on the M/Lj^ ratio and the 
interstellar extinction on the colors. In both these models, the 
burst is of 0.3 Gyr duration, which started around 0.8 Gyr ago. 
Star formation activity in the disk has ceased since then. The 
dotted paths indicate the locus of this point after corrections for 
gas content (fraction of gas mass is 15%) and extinction (length 
of the arrow corresponds to A v = 1.2 mag) are applied. Model B 
successfully reproduces all the observed trends. 



after the burst has stopped, we are forced to truncate the 
lower end of the IMF at 0.4 M Q . Moreover, in order to 
reproduce the observed enhancement with SNIa progeni- 
tors that originate from mass accretion onto white dwarfs 
in binary systems (as it is adopted in our chemical model) 
we need to decrease the slope of the IMF of the massive 
stars to the Salpeter value. Thus, by combining age, lu- 
minosity and dynamical mass of the stellar populations 
we conclude that the disk of M 82 formed very recently 
and with an IMF peaked toward the massive stars. The 
increase in the low-mass cutoff and decrease in the slope 
may be due to a turbulence-enhanced ISM, probably as 
a result of interaction or/and the starburst episode itself. 



5. EVOLUTIONARY STATUS OF M 82 

Disks of normal spiral galaxies are characterized by 
star formatio n that is continuing un -interrupted over the 
Hubble time (|Kennicutt et alJl!994l ). M82 differs in two 
ways from these normal disks: (1) luminosity weighted 
age of the disk is only 0.5 Gyr, which is an order of 
magnitude younger, and (2) there is no present star for- 



mation in the disk. M82 also differs from the normal 
galaxies in the relations between global parameters, it 
is too bright and also met al rich for its low mass value 
ijKarachentsev et alJ l2004j) . The origin of these pecu- 
liarities of M 82 could be understood with our inferred 
SFH, in which almost the entire observed stellar mass 
was formed in a short-duration burst, rather than in 
a continuous manner. In our best-fit model, the disk- 
wide burst started around 0.8 Gyr ago, which is con- 
sistent with the interaction scenario proposed by Sofuc 
(1998). It is interesting to note that a number of in- 
teracting galaxies also show disk-wide star formation 
at present: examples being NGC 4038/9, Arp29 9 and 
NGC 4676 fflibbard & Yunl [19991 iMirabel et all IT99& 
IXu et al] l2000). This kind of widespread star formation 
in galactic disks is probably induced by shocks gener- 
ated in the int erstellar medium following the interaction 
iBarnes Il200l 

Mihos ( 1999) discussed pre-collisional galaxy proper- 
ties that favors star formation in the disk, before the 
burst is turned on in the center. He found that the pre- 
collisional galaxy should be a normal high-surface bright- 
ness galaxy with a big bulge or a gas-rich low surface 
brightness galaxy (LSB). If the pre-collisional M 82 had 
a big bulge, it would have survived the interaction (Sofuc 
1998) . However, the observed bulge of M 82 is small, con- 
sequently, pre-collisional M82 should have been a gas- 
rich LSB galaxy. This is in agreement with our proposed 
SFH. It could be that either M82 was a gas-rich LSB 
galaxy before the interaction, or that it acquired vast 
amounts of gas during the interaction. Below we dis- 
cuss three scenarios that are consistent with the present 
appearance of M82. 

5.1. Kinematically disturbed normal galaxy 

The most important distinguishing characteristic of 
M 82 is its low value of M/L ratio as compared to that of 
normal galaxies. Is this because of an underestimation of 
its mass? The kinematics of the galaxy could have been 
disturbed because of interaction resulting in possible un- 
derestimation of its mass. Blue-band spectra would look 
identical even if we add upto a factor of 3-5 more mass 
in stars older than 5 Gyr as compared to the mass in 
stars 10 times younger. However, the presence of such 
an amount of old stars would have contributed to the en- 
richment of Fe, without the corresponding enrichment of 
a elements, thus resulting in low values of [a/Fe]. Thus 
the observed enhancement of a elements indicates that 
the galaxy really lacks a significant population of stars 
older than 1 Gyr. 

5.2. Rejuvenated star formation in a stripped stellar 

disk 

In the scenario proposed bv lSofuel l)1998f) . much of the 
disk and halo of M 82 were stripped off during its inter- 
action with the members of the M 81 group. However, 
none of the observational studies so far have detected 
the stripped dark matter halo, containing the old disk. 
It would be interesting to carry out a search for such 
a system in the M 81 group. Even in this scenario, the 
gas and stars in the central kiloparsec region are not 
stripped, and hence the observed a enhancement of the 
nuclear gas requires vast amounts of star formation, and 
inflow of the metal-enriched material, after the original 
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disk has been lost. In other words, the SFH we have in- 
ferred is consistent with this scenario, and hence M 82 is 
in the process of re-building its stellar disk. 

5.3. A nearby galaxy in formation? 

Analysis of galaxies at redshifts up to 1 support a mor- 
phological evolution over the past few gigayears. The 
catalyst that drives the morphological evolution is the 
major mergers between galaxies. During the merging 
process galaxies grow both in size and luminosity. Data 
suggest that at least 50% of t he nearby disk galax ies are 
formed through this scenario ijHammer et alJl2005|) . The 
most important ensemble of galaxies that are thought 
to be caught in the formation process are the luminous 
blue compact galaxies (LBCGs) seen between redshifts 
0.2-1.0. The characteristic signature of these galaxies is 
the enhanced st ar formation in their nuclear regions. In 
a recent study, iNoeske et alJ 1)200611 have characterized 
the properties of the stellar disk surrounding the com- 
pact nucleus, and found that they are, in general, short 
in scalelength as compared to the scalelengths of local 
galaxies. 

M 82 resembles these LBCGs in having high luminos- 
ity, short scalelength, a nuclear starburst and low mass. 
Moreover, interaction is responsible for all the observed 
characteristics in M82. Thus in M82, we may be wit- 
nessing stars being formed for the first time in its disk, 
probably in the debris left behind in the interaction. In 
other words, M 82 could be a nearby galaxy in formation. 

6. SUMMARY 

We used spcctrophotomctric, photometric and dynam- 
ical data of the inner «3' (3 kpc) radius of M 82, to study 
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ration of the burst is not enough to reproduce both the 
enhancement of the a elements and the low value of the 
M/Lk ratio, and a suitable top-heavy IMF must be in- 
voked. It is possible that the star formation started in 
the whole disk about 0.8 Gyr ago, induced by shocks 
in the interstellar medium following the interaction with 
members of M 81 group, and now continues only in the 
central regions. 
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